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We demonstrate levitation and three-dimensionally stable trapping of a wide variety of parti¬ 
cles in a vacuum chamber through the use of the thermophoretic force in the presence of a strong 
temperature gradient. Typical sizes of the trapped particles are between 10 pm and 1 mm at a 
pressure between 1 and 10 Torr. The trapping stability is provided by the geometry of the tem¬ 
perature field, as well as the transition between the free molecule and hydrodynamic regimes of the 
thermophoretic force. To quantitatively measure the thermophoretic force, we examine the levita¬ 
tion heights of spherical polyethylene spheres under various experimental conditions and determine 
the temperature gradient needed to levitate the particles. A good agreement between our experi¬ 
mental observations and theoretical calculations is obtained. Our system offers a new platform to 
study thermophoretic phenomena and to simulate dynamics of interacting many-body systems in a 
microgravity environment. 

PACS numbers: 37.10.Mn, 37.10.Pq 


Levitation of particles in ground-based experiments 
provides an ideal platform for the study of their dynam¬ 
ics and interactions in a pristine isolated environment, 
and has broad applications in fundamental and applied 
sciences. Electromagnetically levitated atoms, ions, and 
molecules in vacuum have opened the field of quantum 
degenerate gases [1, 2], cold chemistry [3], and atom- 
based quantum information [4]. Optical levitation of 
micron-size objects [5] provides a new playground to in¬ 
vestigate aerosols in atmosphere [6] and Brownian motion 
of microspheres [7]. 

Levitation of macroscopic particles in low pressure is of 
particular interest due to its wide applications in space, 
atmospheric, and even astro-chemical research. Several 
schemes to levitate glass spheres and water droplets have 
been realized based on optical tweezers [8], high mag¬ 
netic field gradients [9, 10], and electrodynamic balance 
[11]. Experiments on dusty plasma further show many- 
body dynamics of levitating dust particles. An interest¬ 
ing scheme was demonstrated recently to levitate generic 
particles above a hot surface by 100 ^m based on the 
Knudsen compressor effect [12, 13]. The authors further 
demonstrated levitation of ice particles 1 cm above a hot 
surface for a few seconds [14]. 

In this paper, we report an experimental setup to 
achieve stable, three-dimensional trapping of various par¬ 
ticles at pressures of F’=l~10 Torr. The levitation and 
trapping do not require external electromagnetic forces 
or the Knudsen compressor effect, and are achieved with 
only the thermophoretic force in the presence of a strong 
temperature gradient. We observe stable trapping of var¬ 
ious objects including ice particles, ceramic spheres, hol¬ 
low glass bubbles, polyethylene spheres, and lint. The 
stability of levitation allows us to measure the levitation 
heights of the particles as well as the dynamics. Based on 
a numerical model to describe the temperature distribu¬ 
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FIG. 1. Experimental apparatus. The copper plate and liq¬ 
uid nitrogen bucket establish a temperature gradient. The 
pressure is regulated by a vacuum pump and a needle valve 
that opens to atmosphere. Temperature regulation of the 
copper plate is achieved with a Peltier cooler. An absorption 
imaging setup using collimated blue LED light captures the 
shape of levitating particles. A second camera records their 
dynamics. During the levitation, we continuously monitor the 
temperatures of the top and bottom plates. The plates and 
the chamber are electrically grounded to avoid electrostatic 
forces on the levitating particles. 


tion in our system, we obtain a good agreement between 
our experimental observation and theoretical calculation 
over a large range of temperatures where stable levitation 
is achieved. 

Our experimental setup is displayed in Fig. 1. At the 
center of the vacuum chamber is a 10-mm gap between a 
warm copper plate and a stainless steel bucket filled with 
liquid nitrogen at 77 K. The bottom copper plate has a 
diameter of 33.7 mm and is attached to a Peltier cooler 
for temperature control. A thermistor is embedded in the 
plate to allow for temperature monitoring. The stainless 
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steel bucket has a smaller outer diameter of 25.4 mm. 
The bottom surface of the bucket acts as the top plate 
and its temperature is monitored by a K-type thermo¬ 
couple. In the presence of a large temperature gradient, 
particles are levitated and trapped near the mid-plane 
between the the top and bottom plates. 

The levitation comes from the thermophoretic force 
in our chamber. As illustrated in Fig. 2, air molecules 
coming from the hotter side transfer on average a larger 
momentum to the particle than those from the colder 
side. A net momentum transfer due to collisions with 
molecules overcomes gravity and pulls the particles to¬ 
ward the symmetry axis of the two plates. 

In the presence of a temperature gradient VT, the ther¬ 
mophoretic force acting on a sphere of radius a is given 
by [15], 

Fth = -/r —VT, (1) 

V 

where k is the thermal conductivity of the medium (rar¬ 
efied air in our case), v = ^J2k'QT/m is the most proba¬ 
ble speed of the air molecules, m is the molecular mass, 
fee is the Boltzmann constant, and /t is a dimensionless 
thermophoretic parameter that depends on the Knudsen 
number Kn= i/a. The Knudsen number is defined as the 
ratio between the molecular mean free path i = fj,v/P 
and a, where fi is the dynamic viscosity. In the free 
molecule regime where Kn3> 1, /t reaches the maxi¬ 
mum value of 16y^/15 « 1.89. In the hydrodynamic 
limit where Kn<C 1, however, fx approaches zero and 
the thermophoretic force vanishes [15]. 

The thermophoretic force described in Eq. (1) is as¬ 
sociated with heat transfer in the thermal conduction 
regime. In a temperature field T(x), heat flows accord¬ 
ing to Fourier’s law as Q = —nVT. Associated with the 
heat flow, momentum flow per unit area can be estimated 
as p/A = Q/v, where A is the effective cross section of 
a particle. The force on a sphere with radius a is thus 
F = p = —7r(a^K/u)Vr, and is consistent with Eq. (1) 
in the molecular flow regime. 

To stably levitate and trap a particle, one requires 
a restoring force toward the trap center for small dis¬ 
placements in all directions. It may appear that a three- 
dimensional minimum in the temperature field T{x) is 
required. Such condition, however, cannot be satisfied 
without a heat sink under the condition that V • Q = 0 
(Earnshaw’s theorem). Together with Fourier’s law, the 
temperature field in free space satisfies 

V • (kVT) = 0. (2) 

Since the thermal conductivity of air increases with 
temperature [16], the above equation gives V^T = 
— \'\/T\'^k'{T)/k{T) < 0, which excludes the possibility 
of stably confining a particle in free space based solely 
on a temperature gradient. 
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FIG. 2. Simulation of temperature distribution in the cham¬ 
ber is obtained by solving Eq. (2). Boundary conditions hold 
the top plate at 77 K, bottom plate at 277 K, and chamber 
wall at 298 K. Temperature contours are labeled. In the pres¬ 
ence of a temperature gradient, collisions with background 
gas particles exert a net force opposite the direction of the 
temperature gradient owing to the larger momentum transfer 
from hotter molecules. Due to our geometry, particles expe¬ 
rience both an upward force counteracting gravity mg and a 
radial confining force. 

The requirement for stable thermophoretic levitation 
is given by V • + Fc) < 0, where the Fq is the 

gravitational force. Using Eqs. (1), (2), and V • Fiy = 0 
in free space, we obtain the necessary condition for stable 
levitation: 

-V ■Fth = Ka^\VT\‘^4k—>^- (3) 

dl V 

Since thermal conductivity k is positive, the above con¬ 
dition is equivalent to d{fT/VT)/dT > 0. 

An example of the temperature distribution in our sys¬ 
tem is shown in Fig. 2 and is numerically calculated based 
on Eq. (2). Here one can see a positive curvature of tem¬ 
perature in the radial direction that offers radial confine¬ 
ment. In the vertical direction, however, the curvature 
of the temperature is negative, shown in more clearly in 
Fig. 3(a). 

The key to stabilizing particles in the vertical di¬ 
rection is the thermophoretic parameter fx, which is 
strongly suppressed when particles enter the hydrody¬ 
namic regime. Since our system is isobaric p = nksT = 
const., lower temperature near the colder top plate sug¬ 
gests higher density. Particles moving near the top plate 
are thus deeper in the hydrodynamic regime with a cor¬ 
responding smaller thermophoretic force. The total force 
thus switches sign when the particles are over-levitated. 
The distribution of temperature, temperature gradient, 
and levitation force in our system for a typical operation 
are shown in Fig. 3. 
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FIG. 4. Gallery of levitating particles, (a) Low magnification 
image of one polyethylene sphere (Gospheric, UVPMS-BR- 
0.995). High magnification images of (b) one ice particle, (c) 
an aggregate of three hollow glass babbles (Miapoxy, 406354), 
(d) one ceramic sphere (Miapoxy, 406344), (e) one piece of 
lint, and (f) an aggregate of three polyethylene spheres. 


FIG. 3. Thermodynamic qnantities and thermophoretic force 
simnlation for a 10 fim polyethylene sphere with p=998 kg/m^ 
and P=6 Torr. Using results in Fig. 2, the temperature shown 
in (a) and temperature gradient shown in (b) are plotted as 
functions of height, (c) Thermophoretic force (solid line), 
calculated from Eq. (1), is compared with the gravitational 
force on the particle (dashed line). The intersection of the two 
curves indicates the position of the levitating particle. The 
weaker levitation force as the particle rises offers the levitation 
stability. 


The experimental procedure goes as follows. We first 
deposit the desired particles on the bottom copper plate. 
After sealing the chamber, we pump the chamber down 
to F = 1~10 Torr. A pressure gauge (Thyracont, VD8) 
continuously records the pressure throughout the exper¬ 
iment. We then add liquid nitrogen to the steel bucket 
and the top plate quickly reaches 77 K. We launch the 
particles by vibrating the chamber, and once above the 
bottom plate the particles experience the upward ther¬ 
mophoretic force. Depending on the strength of vibra¬ 
tion, the pressure, and the type of particles, the number 
of levitated particles ranges from one to hundreds. Typi¬ 
cally ceramic spheres and hollow glass bubbles tend to be 
levitated in greater numbers than polyethylene spheres. 
Lint is occasionally introduced into the chamber, and we 
observe stably levitating lint as long as 1 mm in length. 
Fig 4 displays a gallery of levitating objects of different 
types. 

The levitation procedure for ice particles differs from 
the above. The top plate is cooled down to 77 K before 
we seal the chamber. Water vapor condensing on the top 
plate creates a layer of ice particles. The chamber is then 
pumped down to the desired pressure. Upon external 
vibration of the chamber ice particles fall from the top 
plate and some are levitated. 


Our detection system consists of two digital cameras 
along with a light-emitting diode (LED) that illuminates 
the levitating particles. One camera with high magni¬ 
fication faces the collimated LED beam and provides a 
closeup view of the particles. Representative images are 
shown in Fig. 4 (b)-(f). The second camera records the 
dynamics of particles with a full field of view including 
the top and bottom plates, shown in Fig. 4 (a). 

To quantitatively characterize the levitating force at 
play, we analyze the measured height of singly levitating 
polyethylene spheres. We focus on single particle levi¬ 
tation to avoid the complex dynamics of multiple parti¬ 
cles. In this case, particles remain stably levitated for 
our longest observation times of one hour, giving us am¬ 
ple time to study the conditions of levitation. To aid 
comparison with theoretical calculations [15], we focus 
on spherical particles. After examination with a micro¬ 
scope, polyethylene spheres were chosen as they appeared 
the closest to ideal spheres. 

Since levitating particles are easily affected by air cur¬ 
rents resulting from pressure changes, we keep pressure 
constant and vary the temperature of the top plate. To 
achieve this we allow the liquid nitrogen in the steel 
bucket to evaporate, which slowly warms up the top plate 
over 5^10 minutes. During this time, we synchronously 
record the temperature of the plates, pressure, and par¬ 
ticle height. 

Examples of singly levitating polyethylene spheres are 
shown in Fig. 5, where the top plate temperature slowly 
increases from 77 K toward room temperature. The par¬ 
ticle’s levitation height gradually decreases as the top 
plate warms up because of the decreasing temperature 
gradient in the chamber. 

Stable levitation is characterized by the balance of 
gravity and the thermophoretic force: 
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FIG. 5. Dynamics of levitating polyethylene spheres dur¬ 
ing a temperature increase, (a) Levitation heights of single 
polyethylene spheres. After levitation we allow the top plate 
to warm up and record the levitation height. Circles (black) 
correspond to a sphere of radius a = 12.7 pm at P = 3.9 Torr, 
triangles (blue) correspond to a particle of a =9.8 pm at P = 
6.4 Torr, and squares (red) correspond to a particle of radius 
9.3 pm at P = 7.3 Torr. (b) Temperatures and temperature 
gradients at the location of particles obtained by simulation 
of temperature distribution in the chamber. Theoretical pre¬ 
dictions from Takata (solid line) [17] and Yamamoto (dashed 
line) [18] are plotted for comparison. 


Fth = pVg, (4) 

where p = 998 kg/m^ is the density of polyethylene 
spheres, g = 9.8 m/s is gravitational acceleration, and V 
is particle volume determined from the particle radius. 
To calculate the thermophoretic force we determine the 
thermodynamic quantities at the location of the parti¬ 
cle, and extract the thermophoretic parameter fx for a 
specific Knudsen number Kn from Eqs. (1) and (3). 

To determine the thermodynamic quantities of air at 
the particle’s location, we solve the air temperature dis¬ 
tribution in the chamber using the steady state heat dif¬ 
fusion equation, Eq. (2), with the appropriate boundary 
conditions. From the distribution we compute the tem¬ 
perature and temperature gradient at the location of the 
particle. To obtain the thermal conductivity n and dy¬ 
namic viscosity p of rarefied air, we obtained their func¬ 
tional forms by fitting the data in Ref. [16] and evalu¬ 
ate them at the position of the particle. We also calcu¬ 
late the Knudsen number Kn= i/a at the particle posi¬ 
tion. Finally, we determine the value of fx for Eq. (3) to 
hold. The procedure is repeated during the entire warm¬ 
ing up process. The measurements on three independent 
polyethylene spheres are shown in Fig. 6. 

Our measurements are consistent with the theories of 
Takata et al. [17] and Yamamoto et al. [18] within ex¬ 
perimental uncertainties. The discrepancies of 5 ^ 10% 


FIG. 6. Comparison between fx from experimental measure¬ 
ments and theoretical models. Thermophoretic parameter fx 
is calculated from the data shown in Fig. 5. The solid line 
represents interpolation of the theoretical calculation of fx 
by Takata et al. [17]. The dashed line represents interpola¬ 
tion of the theoretical calculation of fx by Yamamoto et al. 
[18] with fci 2 = 0.03, where fcr 2 is the ratio of the thermal 
conductivity of gas and the levitating particle. 


between the measurements on three particles are possi¬ 
bly caused by systematic uncertainties of our experiment. 
The particle radii determination suffers from an uncer¬ 
tainty of approximately 1 pm due to our imaging resolu¬ 
tion. Furthermore, the two plates do not have a uniform 
temperature distribution, and the particles can levitate 
slightly away from the symmetry axis. These result in 
uncertainties of temperature and temperature gradient 
at the location of the particle less than 10%. 

Another systematic effect is the photophoretic force 
coming from external radiation. Such force is significant 
in former levitation experiments [14, 19]. We estimate 
the photophoretic force in our system based on Refs. [19, 
20] and find that the radiation power of 10 W/m^ from 
the copper plate at 277 K and 2 W/m^ from the stainless 
steel bucket at 77 K contribute a net levitation force of 
< 0.03 pN on the particles. Compared to the typical 
particle weights of 25~100 pN, the photophoretic force 
is negligible. 

Convective flow is another potential systematic. Con¬ 
vective flow occurs due to the Rayleigh-Benard instabil¬ 
ity, which occurs in a rarified gas at high Rayleigh num¬ 
ber given by [21, 22] 


_ 2048 1 - r 1 
^ “ 757r (1 -b r)2 FrKn*^ 


> 1708. 


( 5 ) 


Here Kn* = i/H is the Knudsen number of the system, 
Fr= /gH is the Froude number, r = T^/Th is the tem¬ 
perature ratio of the two plates, and H is plate separa¬ 
tion. Given typical parameters of our system = 77 K, 
Th = 277 K, i = 7 pm, H = 1 cm, and v = 290 m/s, our 
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Rayleigh number Ra= 9.1 is well below the critical value 
of 1708. We thus conclude that convection is insignifi¬ 
cant. 

In conclusion we demonstrate a versatile experimen¬ 
tal system that levitates particles of various types, and 
stably confines them. We perform a quantitative study 
on singly levitating polyethylene spheres and establish a 
numerical model to determine the thermodynamic quan¬ 
tities at the location of the particles. A direct comparison 
between experiment and theory on the levitation condi¬ 
tion allows us to extract the thermophoretic parameter in 
the range of Knudsen number 0.5 <Kn< 1.0. We report 
good agreement with theoretical models. 

Our system offers a new platform to study levitation 


and interactions of generic particles. When two and more 
particles are levitated, we observe interesting dynamics 
in the presence of interparticle interactions and external 
forces. Our system also offers the possibility to simulate 
dynamics of macroscopic particles in a microgravity envi¬ 
ronment, enhancing our capability to detect feeble forces 
between macroscopic objects in ground-based laborato¬ 
ries. 
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